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Abstract 

' The asymmetries between the spectra of leading and nonleading charmed mesons 

Q\ ■ measured in E~A interactions at pi= 340 GeV/c in the WA89 experiment are 

described in the framework of Quark-Gluon String Model. There are two versions of 
the model under consideration: one of them includes the sea charm quark-antiquark 
■ pairs and the other one does not. It's shown that the asymmetries between D~ and 

^ : d+ -meson spectra and between D s and D^~- meson spectra can be fitted by QGSM 

curves obtained with the same parameters as charm asymmetry in ir~ A experiments 
£NJ ■ described in previous studies. The QGSM results are compared with the calculations 

in the next-to-leading approximation of perturbative QCD approach carried out by 
\ the other authors. 

(N ' 

-xf : 

^ ; 1 Introduction 

OV 

^ ! The asymmetries between the spectra of D and D + , D s and Df mesons were measured 

Ph| recently in H~A interactions at pi= 340 GeV/ c in the WA89 experiment ]l] . It seems to be 

interesting to compare this asymmetries with those ones obtained in the ir~A experiments 
in order to understand the influence of quark composition of beam particle on a heavy 
flavored particle production and to extract some specific features of the strange-charmed 
meson spectra caused by the presence of a strange valence quark in S~ hyperon. The 
difference in xp spectra {xp = 2p n /yfs) of leading and nonleading particles has been 
discussed recently and several theoretical models explained successfully the asymmetry as 
an effect of an interplay between the quark contents of the projectile and of the produced 
hadron. The charmed mesons containing ordinary quarks of the same type as beam 
particle have higher average x value. The asymmetry, defined as: 

At dN D ~/dx-dN D+ /dx 
^ ~ dN D ~ Jdx~+ dN D+ j dx ' ( } 

is a function rising with x. There are two different approaches in the theory for the 
description of this effect. The first one is based on perturbative theory of QCD. It takes 
into account the recombination of 'intrinsic charm' 

quarks with valence quarks of the projectile as an origine of the asymmetry 
||. The other phenomenological models exploit the properties of fragmentation func- 
tions in order to insert the asymmetry . We will not discuss here the details of recombi- 
nation models but we are going to concentrate on a nonperturbative approach known as 
the Quark Gluon String Model (QGSM) Q. This model fits well the leading/nonleading 
charm asymmetry for n~p experiments ||. 



1 



2 The Quark Distributions in QGSM. 

The inclusive production cross section of D-mesons is written as a sum over n-Pomeron 
cylinder diagrams: 



h = x ~^r( s ' x ) = / E ~^ d2p± = £ a ™( s Vn( s > x ) ( 2 ) 



where function y^(s, x ) ls a particle distribution in the configuration of n cut cylinders 
and a n is the probability of this process. The parameter of the supercritical Pomeron used 
here is A P = Op(0) — 1 = 0, 12. The detailed formulae for a n and (p® in pp-interactions 
can be found in || . 

The distribution functions in case of S - -p collisions are: 

tf(a, x) = a^(F^(x + )F^(x.) + F${x + )Fj?\xJ) + 2(n - 1)^1)41 (3) 

where a® is the density parameter of quark- ant iquark chain fragmentation into the given 
type of mesons. 

The particle distribution on each side of chain can be build on the account of quark 



contents of beam particle (x + = {x + \jx 2 + Xj_)/2, x± = 2m^/ sqrts) and of target 
particle (x- = (x — \Jx 2 + x 2 L )/2): 

^ (n) (*+) = \Fl»\x + ) + lF<r\x + ), 

F^(x + ) = + (4) 

F^ n \x-) = ^i n) (x-) + ^FW(x_), 



F%\x_) = \f^(x.) + \fS 
Each F,i(x±) is constructed as the convolution: 



Fi{x ± ) = ^ fU^)—^F(—)dxi, (5) 

J x± X\ X\ 

where f l {x\) is a structure function of i-th quark wich has a fraction of energy x\ 
in the interacting hadron and T>f(z) is a fragmentation function of this quark into the 
considered type of D or D s -mesons. 

The structure functions of quarks in interacting proton have already been described 
in the previous papers 0. In the case of hyperon they depend on the parameter of the 
Regge trajectory of <£>-mesons (ss) because of s-quark contained in £~: 

_ x ^a R {0)-2a N (0)+a R (0)+n-l 
»)- 2 ^(°)(l_ Xl )-^(0)+n-l ; ( 6 ) 
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where 0^(0) =0. The constants C^' are determined by normalization conditions: 



o 



f i {x 1 )dx 1 = 1. 



3 The Fragmentation Functions. 

The fragmentation functions are constructed for the quark and diquark chains according to 
the rules proposed in 0. The following favoured fragmentation function into D~-mesons 
was written for the strange valence quark: 

Vf'(z) = -(1 - 2 )-^(°)+ A (l + a^z 2 ), (8) 

where \=2a' D „ (0)p 2 ±D *. An additional factor (1 + af s z 2 ) provides the parametrization 
of the probability of heavy quark production in the interval z=0 to z — > 1. The values of 
the constant af s will be discussed later. 

The function for the nonleading fragmentation of d-quark chain is: 

V® + (z) = -(1 - z y a R(°)+ x + 2 ( 1 - a R(°))+ A ^ y (9) 

where A^ = or(0) — a^(0). The function of the nonleading fragmentation of the 
diquark chain is the following: 

V® 7 (z) = -(1 - 2)-M°)+*+2(l-aH(0))+A,*. (10) 

where «r(0)=0.5 and = an(0) — aty(0). 

The following fragmentation function corresponds to the version of the diquark frag- 
mentation into D~-mesons: 



V® a (z) = — (l-2) afl(0) ~ 2ajv(0)+A+A ' / '(l + af s z 2 ) + — (i _ /2 ) a «( )~ 2ajv ( ) +A+A v +A * +1 (11) 

4 The Asymmetry Suppression Causes. 

Some fractions of sea quark pairs in hyperon, dd and ss, are to be taken into account 
as far as they suppress the leading/nonleading asymmetry. The structure functions of 
ordinary quark pairs in the quark sea of hyperon can be written by the same way as the 
valence quark distributions: 



f*_( Xl ) = C W x -^(0)(l_ Xl )^(0)-2a JV (0) + A v+ n-l +2 (l- afl ,(0)) ; ( 12) 

where sea quarks and antiquarks have an additional power term 2(1 — or(0)) corre- 
sponding to the quark distribution of two pomeron diagram which includes one sea quark 
pair. 

The structure function for strange sea quarks obeys the same rules: 
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f*_( Xl ) = ^^^-^(^(i-^^^W-a^w+^+n-i^Ci-^co)) (13) 
where A v = Or(0) — a^(0) and 5 StS =0.25 ( see [§]). 

The fragmentation function of strange sea quark( or antiquark) into D s mesons has 
the following form for mesons of the both charges: 

The additional fragmentation parameter a^ s is equal to the fragmentation parameter 
for D-mesons. 

5 The Intrinsic Charm Distribution. 

As soon as we accounted dd and ss fraction in the quark sea of hyperon some fraction 
of charmed sea quark are to be considered as well. This small heavy quark admixture 
plays an important role due to its strong impact on the difference between leading and 
nonleading charmed meson spectra. 

The charmed sea quark structure function is similar to the distribution of strange sea 
quarks: 

/£f(zi) = C '^ ) 5 CiC -xr"' (0) (l-^i) afl(0)_2Qjv(0)+A ^ +n_1+2(1_Qfl(0)) (15) 

where 5 Cj5 is the weight of charm admixture in the quark sea of hyperon. In fact it is 
not neseccarily to be equal to the charmed quark fraction in quark sea of pion ||. This is 
only one parameter we can vary for S~ interaction after the best fit of pion experimental 
data which had been done before. The value of <5 Cjg can be estimated in the description 
of the WA89 data on D s and D meson asymmetries. 

Fragmentation functions are the following: 

V^(z) = -z 1 -^ (1 - z)- aR(0)+x (16) 

and for D s : 

V^(z) = V-^(o)(! _ z )-M0)+a 

6 The Final Plots and The Comparisons. 

The main parameter of QGSM scheme which is responsable for leading/nonleading charm 
asymmetry is a\. It is the parametrization parameter of leading fragmentation function 
dependence on z — > 1. The fraction of charmed sea quarks, <5( Cj c), is the second parameter 
in this calculations which makes the asymmetry lower because of the equal amounts of 
D + and D~ mesons produced by each sea charmed quark pair. Two sets of this couple of 
parameters were chosen in the description of tt~A reaction data: ai=4, 5( Cj c)=0 and a\= 
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Figure 1: Asymmetry between D and D + spectra obtained in E791 (black squares) [0 
and QGSM curves at two versions of parameters: solid line corresponds to the following 
set of parameters a\ = 10,5( Cjg ) = dashed line is a result of QGSM fit with 

a>i = 4, 5( c ,c) = 0. 



10, 5( C) c) = 0.05ya(). We concider here these two values of d\ taking the <5( CjE ) as more or 
less free parameter. 

The two curves displayed in Fig.l represent the fits of E791 pion beam experiment data 
[[J with two sets of parameters discussed above. Data of the WA89 experiment are given 
in Fig. 2 and Fig. 3 with the same parameters. It should be mentioned that the smaller 

fraction of charmed sea quarks was taken into account (5( Cj s) = O.OlyojP for to describe 
both D~ / D + andD~ / asymmetries with the fragmentation parameter a\ = 10. 

The resulting curves obtained in several theoretical models |Tl|, [12| are also shown in 
these Figures as well. 



7 Conclusions. 

There are several conclusions derived from the calculations discussed in the article: 

1) Data of the WA89 experiment on charm production asymmetry can be described 
within the framework of Quark-Gluon String Model with the same asymmetry parameter 
di = 10. as E791 data for tt~A reaction. 

2) D~/D + and D~/D+ asymmetries measured with S~ beam are more sensitive to 
the weight of charmed quark pairs in the quark sea of interacting hyperon (<$( Ci c)=0.01) 
than it could be seen at ir~ beam interaction (5( cg )=0.05). 

3) D~ / asymmetry is higher than D~ / D + asymmetry because strange quark pairs 
suppressing the asymmetry at D s production have lower weight in quark sea of hyperon 
than ordinary dd pairs which cause the suppression of D~ / D + meson asymmetry. 

4) The both of charmed meson asymmetries are nonzero values at xp=0 in these 
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Figure 2: D / D + asymmetry measured in WA89 and theoretical calculations: solid line 
corresponds to the following set of QGSM parameters a\ = 10, <5(C, c) = 0.01^0^; dashed 
line is a result of QGSM fit with a\ = 4, 5i C £) 



0; dashed-dotted line is the result of ITT 



and dotted line corresponds to A(x) predicted in [jl2 



calculations at WA89 energy and diminish with the increasing of energy. 
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Figure 3: D s / asymmetry measured in WA89 [p]]; theoretical curves are the same as 
in Fig. 2. 
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